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Abstract 

Aluminum-metal  oxide  energetic  compositions  with  components  mixed  on  the  nano-scale  are  substantially  more  reactive 
than  conventional  thermites  and  are  of  interest  as  potential  additives  to  propellants,  explosives,  and  pyrotechnics.  For 
such  nanocomposite  materials  prepared  by  Arrested  Reactive  Milling,  the  exothermic  reactions  leading  to  ignition  were 
detected  to  begin  at  relatively  low  temperatures.  These  materials  are  prepared  by  mechanical  processing  at  room 
temperature,  and  the  nature  of  the  interface  present  between  aluminum  and  the  oxidizer  (metal  oxide,  e.g.,  CuO,  M0O3, 
BFO3,  etc.)  is  unknown.  Initial  estimates  show  that  the  reaction  rates  observed  in  thermo-analytical  experiments  cannot 
be  explained  assuming  that  the  metal  and  oxidizer  are  separated  by  a  substantial  oxide  layer,  similar  to  that  present  on 
the  surface  of  conventional  aluminum  particles.  Experiments  using  a  Thermal  Activity  Monitor  (TAM  III)  quantify  the 
reaction  rates  between  aluminum  and  oxidizer  at  temperatures  as  low  as  30  °C;  the  results  of  these  measurements  are 
presented  and  discussed  in  the  context  of  the  redox  reactions  leading  to  ignition  in  such  materials. 


Introduction 

Multiple  reactive  compositions  include  aluminum 
powders  or  nanopowders,  e.g.,  [1  -  3].  Aluminum 
particles  can  be  mixed  with  various  oxidizers,  including 
ammonium  perchlorate,  metal  oxides,  and  other 
compounds.  In  all  such  mixtures,  the  initial  reactions 
leading  to  aluminum  ignition  are  rate  limited  by  relatively 
slow  diffusion  of  the  reacting  species  through  the 
protective  natural  aluminum  oxide  layer  always  present 
on  the  A1  surface  exposed  to  air  [4-7].  However,  the 
reactive  interfaces  existing  in  nanocomposite  materials 
prepared  by  Arrested  Reactive  Milling  (ARM)  [8]  are  not 
exposed  to  gaseous  oxidizer  during  the  material 
preparation.  Indeed,  during  the  mechanical  milling 
performed  under  inert  environment,  when  a  fresh 
aluminum  surface  is  produced,  it  is  immediately  pressed 
against  that  of  the  solid  oxidizer  particles,  e.g.,  CuO, 
M0O3,  etc.  As  a  result,  fully  dense,  nanocomposite 
structures  are  formed,  as  that  shown  in  Fig.  1.  It  is 
unclear  whether  a  passivating  AFO  ;  layer,  similar  to  that 
present  on  any  aluminum  surface  exposed  to  air,  also 
exists  between  inclusions  or  layers  of  metallic  A1  and 
solid  oxidizer  fully  consolidated  within  each  powder 
particle  prepared  by  ARM.  If  such  an  oxide  layer  is 
present,  the  kinetics  of  the  initial  exothermic  reactions 
leading  to  ignition  are  expected  to  be  similar  to  those  in 
the  composite  materials  comprising  blends  of  aluminum 
powders  with  respective  oxidizers.  However,  if  the 
alumina  layer  is  not  formed,  or  if  the  structure  and 
morphology  of  the  layer  separating  metallic  aluminum 
and  oxidizer  are  different  from  those  of  amorphous 
alumina,  it  is  reasonable  to  expect  that  the  ignition 
kinetics  for  ARM-prepared  materials  are  qualitatively 
different  from  those  of  the  mixed  powders  or 
nanopowders  with  the  same  bulk  composition.  The 


objective  of  this  effort  is  to  establish  whether  initial 
exothermic  reactions  in  the  ARM-prepared 
nanocomposite  thermites  can  be  described  using  a 
conventional,  thermally  activated  diffusion  or  whether  a 
different  reaction  mechanism  should  be  introduced. 


Onset  temperatures  for  the  exothermic  reactions  in 
thermite  compositions 

Multiple  literature  reports  are  available  on  onset 
temperatures  and  reaction  kinetics  for  exothermic 
reactions  in  thermite-type  composite  and  nanocomposite 
materials  prepared  by  mixing  of  nanopowders  and  by 
layered  deposition  of  nanofoils.  The  deposition  is 
performed  at  elevated  temperatures  that  support  a  more 
extended  reaction  between  the  deposited  layers. 
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For  variety  of  aluminum  based  thermites,  including  Al- 
M0O3  [9],  Al-Fe203  [10],  AI-FO5  [11]  and  other  systems, 
the  onset  of  exothermic  redox  reactions  was  detected 
between  500  and  600  °C,  with  lower  temperatures 
observed  for  finer  powders  and  nanopowders.  Similarly, 
the  onset  of  exothermic  reactions  in  the  same  temperature 
range  was  observed  for  Al-CuO  nanofoils  [12,  13].  Note 
that  this  temperature  range  effectively  coincides  with  the 
temperatures  at  which  the  first  oxidation  step  for 
aluminum  powders  and  nanopowders  is  observed  in 
gaseous  oxidizers  [4-7]. 

Detailed  studies  of  exothermic  reactions  in  the  ARM- 
prepared  nanocomposite  thermites  [14  -  17]  showed  that 
the  reaction  onset  is  observed  at  a  much  lower 
temperatures,  between  150  and  250  °C.  This  significant 
difference  in  the  onset  temperature  for  the  exothermic 
reaction  suggests  that  the  reaction  mechanism  and 
respective  rate-limiting  processes  might  be  different  for 
ARM-prepared  materials. 

Focusing  on  low  temperature  processes,  in  this  work 
the  heat  release  during  reaction  in  ARM  material 
2Al+3CuO  in  isothermal  conditions  was  studied  using  a 
microcalorimeter  TAM  111  by  TA  Instruments. 
Experiments  with  2  Al+3CuO  in  TAM  were  conducted  at 
three  different  temperatures:  303  K,  323  K,  and  373  K. 
The  samples  used  in  the  experiments  differed  by  the  aging 
time  and  conditions  prior  to  the  experiments.  All 
samples  were  prepared  by  milling  in  Ar  atmosphere  using 
hexane  as  a  process  control  agent  (the  general  procedure 
is  described  in  detail  elsewhere  [14]),  and  dried 
immediately  after  preparation  for  about  18  hrs  in  Ar 
atmosphere.  After  that,  the  storage  conditions  varied  for 
different  samples.  Table  1  summarizes  the  conditions  at 
which  the  material  was  stored  before  microcalorimetry 
measurements  using  TAM  III. 

Figure  2  shows  the  results  of  TAM  III  measurements 
in  terms  of  both  integral  heat  release  and  heat  release  rate, 
normalized  by  sample  mass.  All  samples  were  exposed  to 
the  furnace  temperature  at  the  time  equal  to  0  s;  however, 
the  measurements  are  unavailable  for  the  initial  period  of 
time  while  the  sample  is  thermally  equilibrating.  The 
traces  are  shown  after  the  TAM  111  signal  became 
meaningful.  The  curves  are  grouped  by  the  temperature 
at  which  the  TAM  III  measurements  were  performed. 
Each  curve  within  a  group  belonging  to  the  same 
temperature  corresponds  to  the  sample  with  different  pre¬ 
treatment  history.  Note  that  only  one  curve  is  shown  for 
experiment  performed  at  303  K,  corresponding  to  the 
freshly  prepared  sample.  For  that  sample,  a  noticeable 
heat  release  is  measured  even  at  this  relatively  low 
temperature.  Samples  stored  for  a  year  after  preparation 
did  not  exhibit  any  reaction  when  treated  at  303  K. 

The  curves  are  seen  to  group  based  on  the  processing 
temperature,  indicating  that  the  processing  temperature 
had  the  major  effect  on  the  reaction  rate,  while  the  pre¬ 
heating  history  resulted  in  measurable,  but  smaller  effects. 


Some  difference  between  the  measured  reaction  rates  was 
observed  for  samples  B  and  D  with  different  storage 
conditions  that  were  processed  at  323  K  with  the  sample 
D  stored  in  Ar  showing  a  slightly  greater  heat  release. 
Reaction  rates  were  effectively  identical  for  the  samples  B 
and  C  that  were  stored  at  the  same  conditions,  in  air. 
Taking  into  account  that  these  samples  were  from 
different  batches  of  powder,  this  suggests  good 
reproducibility  of  the  TAM  III  measurements.  Difference 
for  the  curves  of  the  sample  processed  at  373  K  is  due  to 
preprocessing  of  one  of  them  (E)  at  323  K,  and  its 
reaction  at  373  K  is  slower  than  that  of  the  sample  F  kept 
in  air  at  room  temperature. 


Table  1.  Aging/storage  conditions  for  sample  used  in 
TAM  studies 


Sample 

ID 

Temperature, 

K 

Aging  Conditions 

A 

303 

Used  immediately  after 
drying 

B 

323 

Stored  for  1  day  in  air  at 
room  temperature 

C 

323 

Stored  for  1  day  in  air  at 
room  temperature 

D 

323 

Stored  for  1  day  in  Ar  at 
room  temperature 

E 

373 

Stored  for  1  day  in  air  at 
room  temperature;  Exposed 
to  323  K  for  ~2.5  d 

F 

373 

Stored  for  5  days  in  air  at 
room  temperature 

Time,  hours 

0  10  20  30  40  50 


Time,  s 

Fig.  2.  Heat  and  rate  of  heat  of  reaction  of  2Al+3CuO 
ARM  nanocomposites  at  various  temperatures  obtained 
by  microcalorimeter  measurements. 
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Processing  of  the  microcalorimetry  and  differential 
scanning  calorimetry  measurements 

Differential  scanning  calorimetry  (DSC)  and  TAM  III 
data  (microcalorimetry)  provide  quantitative  measures  of 
the  heat  flow  generated  in  the  sample  upon  its  heating 
(DSC)  or  in  isothermal  conditions  (microcalorimetry). 
This  heat  flow  is  used  to  estimate  the  extent  of  the  redox 
reaction  that  should  occur  to  account  for  the  measured 
heat  release.  Alternatively,  the  values  of  heat  flow  can  be 
translated  into  equivalent  oxide  layer  thicknesses  that 
should  be  grown  as  a  result  of  the  redox  reactions.  As 
described  below,  both  measured  heat  flows  and  inferred 
values  for  the  oxide  thickness  can  be  compared  to  those 
predicted  theoretically  considering  different  reaction 
mechanisms.  Using  this  approach,  the  mechanism  of  the 
initial  exothermic  reaction  can  be  elucidated. 

The  estimates  for  the  growing  oxide  layer  thickness 
were  performed  for  the  experimental  data  described  above 
for  microcalorimetry  measurements  and  for  data  reported 
in  ref.  [14]  for  DSC  measurements. 

The  predicted  thickness  of  the  oxide  grown  is  a 
function  of  the  assumed  sample  geometry.  Two 
geometries  are  considered.  In  the  spherical  geometry,  a 
spherical  oxide  inclusion  is  surrounded  by  A1  matrix.  In 
the  planar  geometry,  layers  of  oxide  and  A1  are  placed  on 
top  of  each  other. 

Figure  3  shows  the  oxide  thickness  estimated  from 
the  DSC  measurements  [14]  for  both  geometries  and 
assuming  that  100  nm  is  the  initial  inclusion  diameter  and 
oxide  layer  thickness  for  the  spherical  and  planar 
geometries,  respectively.  Note  that  for  the  spherical 
geometry,  the  result  of  the  estimate  also  depends  on  the 
thickness  of  the  initial  oxide  layer  assumed  to  be  present 
on  surface  of  aluminum  particle.  This  dependence  is 
relatively  as  long  as  the  thicknesses  of  oxide  layers  are 
small  compared  to  the  particle  diameter.  For  consistency, 
several  initial  values  for  the  initial  alumina  thickness  were 
considered  for  both  geometries. 

It  is  interesting  that  a  much  greater  oxide  layer  should 
be  growing  on  the  flat  surface  (about  6  nm)  as  compared 
to  spherical  particles  (only  about  1.5  nm)  to  explain  the 
heat  effect  measured  upon  heating  the  nanocomposite 
sample  in  DSC  furnace  to  600  K. 

Examples  of  similar  estimates  based  on  the 
microcalorimetry  measurements  are  shown  in  Fig.  4.  As 
in  Fig.  3,  the  calculations  are  performed  for  both  spherical 
and  planar  sample  geometries.  In  this  case,  the  thickness 
of  the  initial  alumina  layer  is  kept  constant  (0.5  nm)  and 
the  diameter  of  the  inclusion  or  layer  thickness  are  varied 
from  100  to  1000  nm  for  both  geometries.  As  for  the 
DSC  measurements,  a  much  thicker  alumina  layer  is 
predicted  to  grow  for  the  flat  geometry  than  for  spherical 
inclusions  to  explain  the  measured  heat  release. 


Formulation  of  a  theoretical  reaction  model 

The  results  of  TAM  III  measurements  were  interpreted 
using  two  reaction  models.  One  of  the  models  was 
conventional  thermally  activated  diffusion  [18],  and  the 
other  was  so-called  Mott-Cabrera  model  developed  to 
describe  growth  of  initial  thin  oxide  layers  at  low 
temperatures  [19,  20].  In  the  latter  model,  the  mass 
transfer  rates  are  accelerated  by  electrostatic  interaction  of 
the  diffusing  ions  and  static  electric  charges  developed 
principally  due  to  the  difference  in  mobility  between  ions 
and  electrons.  Because  the  focus  of  the  present  effort  was 
to  describe  the  onset  of  the  heterogeneous  exothermic 
reactions,  the  oxidation  results  in  the  growth  of  the  initial, 
amorphous  alumina  layer.  The  polymorphic  phase 
changes  in  alumina  occur  at  higher  temperatures  and/or 
substantial  oxide  thicknesses  [5  -  7]  and  thus  play  only 
secondary  role  for  the  processes  of  interest  in  this  study. 
Transport  (diffusion)  of  oxygen  from  the  oxide  core  to 
aluminum  matrix  through  a  growing  layer  of  alumina  was 
considered  as  the  reaction  rate-limiting  process.  Note 
that  this  analysis  is  effectively  equivalent  to  a  case  of 
aluminum  ions  diffusing  to  oxidizer;  in  the  present 
simplified  analysis,  the  nature  of  the  diffusing  species  is 
not  affecting  the  model  outcome. 

Because  reactions  of  reduction  of  the  metal  oxidizers 
used  in  thermite  systems  may  not  occur  in  one  step,  each 
reaction  model  should  represent  a  specific  oxidizer.  No 
change  in  the  general  model  layout  is  required  to  modify 
it  for  a  different  oxidizer;  however,  specific  reaction 
steps,  products,  and  subproducts  should  be  altered.  As  a 
first  step  the  model  formulated  in  this  study  considers  an 
Al-CuO  thermite.  This  material  system  is  selected 
because  of  two  primary  reasons: 

a.  There  is  only  one  known  intermediate  copper 
oxide,  Cu20  that  can  be  formed  upon  initial 
reduction  of  CuO  before  metallic  copper  is 
produced.  For  most  of  the  other  oxides  used  in 
practically  interesting  thermite  compositions, 
more  intermediate  oxides  can  be  formed 
resulting  in  a  greater  number  of  the  anticipated 
reaction  steps. 

b.  The  cross-sectioned  Al-CuO  composites  show 
good  contrast  in  the  SEM  images  and  thus  a 
relatively  straightforward  measurement  of  the 
sizes  of  A1  and  CuO  domains  present  in  the 
actual  samples  is  possible,  as  needed  for 
calibration  and  validation  of  the  proposed  model. 

Schematically,  the  reactant  configurations  used  in  the 
model  are  illustrated  in  Fig.  3.  For  a  spherical  geometry, 
an  inclusion  of  CuO  is  considered  to  be  embedded  into  an 
A1  matrix.  For  a  planar  geometry,  flat  layers  or  CuO  and 
A1  are  placed  on  top  of  each  other.  For  all  cases,  the 
reacting  components,  A1  and  CuO,  are  assumed  to  be 
initially  separated  by  an  amorphous  alumina  layer,  similar 
to  that  present  on  surface  of  any  aluminum  particle 
exposed  to  air.  This  oxide  layer  can  transform  into 
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different  alumina  polymorphs  at  elevated  temperatures 
and  increased  thickness  [5-7]. 

It  is  assumed  that  the  rate  of  oxygen  diffusion  within 
the  oxide  core  is  significantly  higher  than  the  rate  of 
oxygen  (or  aluminum)  diffusion  through  a  growing  layer 
of  alumina.  Therefore,  the  concentration  of  oxygen 
inside  the  oxide  inclusion  is  always  spatially  uniform.  In 
general  case,  the  model  considers  two  possible 
mechanisms:  reduction  steps  for  CuO  through  formation 
of  Cu20  and  Cu,  and  reduction  of  CuO  directly  to  Cu. 
According  to  the  Cu-O  phase  diagram  [21],  at  low 
temperatures  (<  375  °C)  used  in  the  microcalorimetry 
experiments,  only  CuO  and  Cu  are  stable  in  the  binary 
Cu-0  system. 

In  the  case  of  conventional  thermally  activated 
diffusion,  the  growth  of  the  alumina  layer  thickness,  /?,  is 
described  by  the  ‘parabolic’  law  in  case  of  planar 
geometry: 

d±  =  D^  (1) 

dt  °  h 

where  dh/dt  is  the  rate  of  growth  of  the  layer  thickness, 

CQ_  is  the  oxygen  ion  concentration  in  oxide  inclusion; 

and  oo  is  the  diffusion  coefficient  for  oxygen  ions.  For 


the  spherical  geometry,  the  equation  takes  form: 


dt  0 


C 


1 


(2) 


Al20, 


where  VCuO  and  rAi2Q3  are  the  radii  of  the  oxide  inclusion 


and  alumina  layer,  respectively.  The  radii  are  changing  as 
the  oxidation  continues,  considering  the  amounts  of  the 
reacting  species  and  respective  densities  of  Al,  CuO/Cu20 
mixture,  and  A1203. 

For  the  thermally  activated  diffusion,  diffusion 
coefficient  is  described  by  an  Arrhenius  equation: 


Da_  =  D  ■  exp 


(3) 


where  R  is  the  universal  gas  constant,  T  is  temperature, 
and  the  pre-exponent  D  and  activation  energy,  E,  can  be 
taken  from  the  aluminum  oxidation  studies  or  treated  as 
adjustable  parameters  of  the  model. 

The  concentration  of  the  diffusing  species  (oxygen 
anions)  at  the  interface  of  aluminum  and  the  growing 
alumina  layer  was  assumed  to  be  zero.  Two  different 
approaches  were  considered  to  determine  the  oxygen 
concentration  inside  the  oxide  inclusion,  C0_.  In  one 
approach,  the  concentration  of  the  mobile  oxygen  anions 
was  calculated  as  equal  to  the  oxygen  concentration  in  the 
gas  in  equilibrium  with  the  mixed  CuO/Cu20  oxides  at  a 
specific  temperature.  These  concentrations  are  calculated 
using  the  p-T-x  Cu-0  phase  diagram  [21].  This  approach 
produces  a  very  small  oxygen  concentration,  resulting  in 
very  low  overall  diffusion  rates.  Another  approach 
assumes  that  the  entire  oxygen  concentration  present  in 


the  CuO/Cu20  oxide  mixture  is  available  for  diffusion 
and  produces  an  upper  limit  for  the  possible  concentration 
of  oxygen  anions.  In  this  latter  case,  the  calculation 
begins  with  the  concentration  of  oxygen  in  pure  CuO  and 
continues  to  adjust  the  oxygen  concentration  as  the  redox 
reaction  removes  oxygen  from  the  copper  oxides  and 
transfers  it  to  alumina. 

An  alternative  reaction  description  for  aluminum 
oxidation  is  given  by  Mott-Cabrera  model  [19,  20], 
considering  formation  of  the  initial,  very  thin,  oxide 
coverage  on  the  fresh  aluminum  surface.  In  this  model, 
the  diffusion  is  substantially  assisted  by  electric  fields 
formed  due  to  the  difference  in  the  mobility  of  oxygen 
ions  and  electrons  from  the  metal  through  the  growing 
oxide  layer.  Electrons  from  metal  diffuse  through  the 
oxide  layer  to  the  oxidizer  side  and  assist  in  the  formation 
of  oxygen  anions.  Production  of  such  anions  assists  in 
diffusion  of  metal  cations  through  the  growing  oxide  and 
their  reaction  with  the  oxygen  anions.  This  reaction  is 
much  faster  than  conventional  diffusion  and  results  in  a 
rapid  growth  of  an  initial  oxide  layer.  The  reaction  slows 
down  at  a  specific  limiting  thickness  of  the  oxide  film, 
with  further  oxidation  governed  by  regular  diffusion.  The 
Mott-Cabrera  electric  field-assisted  oxidation  model  leads 
to  the  following  equation  for  the  growth  rate  of  the  oxide 
thickness  h  [19,  20]: 

A^expf^'H  (4) 

dt  ^  kBT  y 

where  kB  is  the  Boltzmann  constant,  U  is  the  rate  - 
limiting  energy  barrier,  q  is  the  charge  of  migrating 
cations,  a  is  the  distance  between  the  positions  of  the 
potential  well  and  barrier,  Q  is  the  volume  of  oxide 
formed  per  one  cation,  V  is  the  attempt  frequency  of  the 
cation  jump,  n  is  the  number  of  cations  per  unit  area  that 
may  jump  through  the  energy  barrier  U.  EM.C  is  the 
electric  field  generated  in  the  oxide  film  that  assists 
aluminum  cations  to  overcome  the  potential  barrier.  The 
electric  field  for  the  plane  geometry  is  given  by: 

E  =  -V-  (5) 

h 


For  the  spherical  geometry,  electric  field  is  calculated  as: 

£  _  Yi rC"Q  (6) 

rA,  n  h 


The  chemical  reactions  considered  to  evaluate  the 
reaction  enthalpy  compare  the  theoretical  models  with  the 
DSC  and  microcalorimetry  experiments  are 

3CuO  +  2AI  — »  3Cm  +  Al203  (7) 

at  temperatures  <  375  °C,  and 
6CuO  +  2AI  — >  3Cu20  +  AI203  (8) 

at  higher  temperatures. 
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Fig.  3.  Thickness  of  the  alumina  oxide  layer  predicted  to  grow  on  the  surface  of  a  100-nm  CuO  inclusion  in  DSC 
experiments  [14]  considering  a  spherical  inclusion  (left),  and  a  flat  layer  (right). 
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Fig.  4.  Thickness  of  the  alumina  oxide  layer  estimated  to  grow  in  TAM  III  experiments  at  different  temperatures.  The 

calculations  are  performed  for  spherical  inclusions  (left)  and  flat  layers  (right)  with  the  dimensions  varied  from  100  to 

1000  mu.  The  initial  oxide  layer  is  assumed  to  be  0.5  nm  for  all  cases. 
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Model  predictions  and  their  implications 

Specific  predictions  of  the  oxidation  rates  based  on  a 
conventional  thermally  activated  diffusion  model  used  the 
reaction  parameters  identified  in  ref.  [7]  for  the  first 
detected  step  of  oxidation  of  aluminum  powders  in 
thermogravimetric  experiments.  In  ref.  [5  -  7],  the 
oxidation  was  actually  detected  at  higher  temperatures 
than  in  the  experiments  involving  reactive  nanocomposite 
materials  (refs.  [14  -  17]  and  present  TAM  III 
measurements),  so  that  the  Arrhenius  term  introduced  in 
ref.  [7]  was  extrapolated  to  describe  reaction  at  lower 
temperature.  As  mentioned  above,  the  concentration  of 
oxygen  available  for  diffusion  was  taken  to  be  equal  to 
that  in  the  condensed  CuO.  This  oxygen  concentration  is 
substantially  higher  than  that  available  in  experiments 
described  in  refs.  [5  -  7]  resulting  in  a  significant  increase 
in  the  oxidation  rate. 

For  processing  DSC  experiments  presented  in  ref. 
[14],  the  measured  heat  flows  were  directly  compared  to 
those  predicted  using  the  conventional  thermally  activated 
diffusion  model.  Results  of  this  comparison  are  shown  in 
Fig.  6.  The  bottom  plots  show  measured  and  calculated 
heat  flows  and  the  top  plots  show  oxide  thickness 
estimated  to  be  required  to  describe  the  experimental  data. 
Note  that  separate  thicknesses  for  the  initially  growing 
amorphous  oxide  and  for  the  y-Al203  polymorph  formed 
later  on  are  shown,  as  calculated  in  the  oxidation  model. 

The  shape  of  the  experimental  DSC  curves  could  not 
be  reproduced  despite  treating  both  activation  energy  and 
pre-exponent  in  the  oxidation  model  as  adjustable 
parameters.  Even  assuming  the  maximum  possible 
concentration  of  oxygen  inside  the  CuO  inclusion  or 
assuming  an  unrealistically  high  pre-exponential  factor 
and/or  unrealistically  low  activation  energy  resulted  in  an 
under-predicted  reaction  rate  at  low  temperatures,  at 
which  the  onset  of  the  exothermic  reaction  was  detected 
experimentally.  When  adjustable  parameters  were 
selected  in  unrealistic  ranges,  to  predict  at  least  some 
reaction  at  low  temperatures,  the  predicted  reaction,  while 
starting  late,  occurred  at  the  rates  much  higher  than 
experimentally  observed,  indicating  the  critical  deficiency 
of  the  diffusion-based  oxidation  model. 

It  is  also  interesting  that  the  predicted  evolution  of 
the  oxide  layer  is  qualitatively  different  from  that 
presented  in  Fig.  3,  for  which  the  same  experimental 
results  were  processed  without  considering  any  specific 
oxidation  model. 

In  summary,  it  appears  to  be  impossible  to  match  the 
measured  DSC  traces  using  the  thermally  activated 
diffusion  model. 

For  the  isothermal  heating  (TAM  III  experiment),  the 
results  of  experiments  are  compared  to  the  model 
predictions  in  terms  of  the  rate  of  oxide  growth.  The  rates 
of  oxide  growth  for  the  results  presented  in  Fig.  4  are 
shown  in  Fig.  7  as  solid  lines.  Only  results  for  samples 
A,  B,  and  F  are  presented  to  simplify  the  discussion.  As 


noted  above,  the  temperature  of  the  furnace  affected  the 
rate  of  oxidation  most  substantially,  with  only  minor 
differences  caused  by  the  pre-heating  and  storage  history. 

Long-dashed  lines  in  Fig.  7  show  the  oxide  growth 
rates  predicted  using  a  thermally  activated  diffusion 
model.  These  rates  are  substantially  lower  than 
experimental,  despite  using  the  very  high  concentration  of 
oxygen  ions  in  CuO,  as  described  above.  Because  the 
presented  computed  oxidation  rates  are  very  low, 
effectively  no  measurable  growth  occurs  over  the  time  of 
experiment  and  the  changes  in  the  oxidation  rate  are  too 
small  to  be  noticeable  in  Fig.  7.  Again,  this  behavior  is 
qualitatively  different  from  that  observed  experimentally. 
The  oxide  growth  rates  computed  using  the  Mott-Cabrera 
model  are  shown  in  Fig.  7  as  short-dashed  lines.  All 
parameters  describing  the  model  were  taken  from  ref.  [22] 
and  are  suitable  for  oxidation  of  fresh  A1  surface  in 
gaseous  oxidizer.  It  is  likely  that  these  parameters  need  to 
be  adjusted  for  the  interface  growing  between  aluminum 
and  metal  oxide.  It  is  interesting  that  the  discrepancy 
between  the  predicted  and  experimental  data  is  much  less 
than  for  the  diffusion  model.  Most  importantly,  the  Mott- 
Cabrera  model  predicts  the  change  in  the  oxide  growth 
rate  that  is  qualitatively  similar  to  that  observed  in 
experiments.  It  has  been  observed  that  considering  some 
of  the  oxidation  model  parameters  as  adjustable,  the 
experimental  data  can  be  matched  with  high  accuracy. 
This  was  not  possible  with  any  choice  of  adjustable 
parameters  for  the  simple  diffusion  model.  Therefore,  it 
is  concluded  that  the  oxidation  reaction  in  nanocomposite 
thermites  is  initially  well  described  by  the  Mott-Cabrera 
model,  so  that  the  transport  of  the  reacting  ions  is  assisted 
by  electric  field  formed  across  the  growing  alumina  layer. 

Conclusions 

New  microcalorimetry  experiments  as  well  as  analysis  of 
the  previously  reported  DSC  measurements  show  that  the 
onset  of  exothermic  reaction  in  nanocomposite  materials 
prepared  by  ARM  occurs  at  substantially  lower 
temperatures  than  in  other  Al-based  composite  and  nano¬ 
composite  materials.  The  analysis  of  the  experimental 
data  showed  that  a  simple  model  of  thermally  activated 
diffusion  cannot  adequately  describe  the  onset  of 
oxidation  in  nanocomposite  materials  prepared  by  ARM. 
An  alternative  reaction  description  for  aluminum 
oxidation  is  given  by  Mott-Cabrera  model  [19,  20], 
considering  formation  of  the  initial,  very  thin,  oxide 
coverage  on  the  fresh  aluminum  surface.  In  this  model, 
the  diffusion  is  substantially  assisted  by  electric  fields 
formed  due  to  the  difference  in  the  mobility  of  oxygen 
ions  and  electrons  from  the  metal  through  the  growing 
oxide  layer.  Electrons  from  metal  diffuse  through  the 
oxide  layer  to  the  oxidizer  side  and  assist  in  the  formation 
of  oxygen  anions.  Production  of  such  anions  assists  in 
diffusion  of  metal  cations  through  the  growing  oxide  and 
their  reaction  with  the  oxygen  anions.  This  reaction  is 
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Fig.  5.  Schematic  diagram  of  reaction  geometries  considered  in  the  models:  a  single  spherical  inclusion  of  oxidizer  in  an 
aluminum  matrix  (left)  and  planar  layers  (right). 


Fig.  6.  Bottom  plots:  experimental  DSC  trace  (black)  [14]  compared  to  the  calculated  heat  flow  traces  using  the 
developed  model.  The  top  plots  show  evolution  of  the  amorphous  oxide  layer  and  its  transformation  into  7-AFO3.  The 
left  plots  show  a  calculation  with  realistic  diffusion  coefficient.  The  right  plot  shows  a  calculation  with  unrealistically 
high  pre -exponent  for  the  diffusion  coefficient  selected  to  shift  the  reaction  onset  to  a  lower  temperature. 
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Fig.  7.  Rates  of  A1  oxide  growth  estimated  from  experimental  microcalorimery  data  (solid  lines)  and  their  comparison 
with  the  theoretical  predictions  using  thermally  activated  diffusion  model  (long  dash)  and  Mott-Cabrera  model  (short 
dash)  for  spherical  and  planar  geometries.  Initial  A1203  layer  is  assumed  to  be  0.5  nm  thick  for  all  cases. 
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faster  than  conventional  diffusion  and  results  in  a  rapid 
oxide  growth  at  low  temperatures.  The  reaction  slows 
down  at  a  specific  limiting  thickness  of  the  oxide  film, 
with  further  oxidation  governed  by  regular  diffusion. 
Preliminary  estimates  suggest  that  both  the  experimental 
DSC  curves  and  microcalorimetry  measurements  can 
indeed  be  quantitatively  described  using  the  Mott-Cabrera 
reaction  model. 
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